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Abstract 

We discuss the effects of a possible time variation of the Newton constant Gn on light 
elements production in Big Bang Nucleosyntesis (BBN). We provide analytical estimates 
for the dependence of primordial abundances on the value of the Newton constant during 
BBN. The accuracy of these estimates is then tested by numerical methods. Moreover, 
we determine numerically the response of each element to an arbitrary time-dependent 
modification of the early universe expansion rate. Finally, we determine the bounds 
on possible variations of Gn which can be obtained from the comparison of theoretical 
predictions and observational data. 


‘E-mail: bambi@fe.infn.it 
^E-mail: giannotti@fe.infn.it 
^E-mail: villante@fe.infn.it 


1 




1 Introduction 


The idea that fundamental constants may vary with time dates back to Dirac pQ . Even if, at 
present, there is no robust experimental evidence in favor of this possibility, this idea contin¬ 
ues to be widely discussed, since many extensions of the standard theories (e.g. superstring 
theories, scalar-tensor theories of gravitation, etc.) predict a variation of ’’fundamental 
constants” with time. 

On a pure theoretical level, the space-time dependence of the fundamental parameters^ 
is forbidden by the Strong Equivalence Principle (SEP), and in particular by the statement 
of Local Position Invariance (LPI) (see e.g. [2]). On the other hand, the weakest form 
of the Equivalence Principle, the so-called Einstein Equivalence Principle (EEP) - which 
is the essence of the geometrical theory of gravitation ~ refers (in its LPI) only to non- 
gravitational physics and, thus, allows the Newton constant Gn to be time-dependent. This 
is what happens, e.g., in the Brans-Dicke, or more generally in the (multi-) scalar-tensor 
theories of gravity. 

The above argument suggests that Gn has a special role in the subject of time variation 
of the fundamental parameters. A dependence of Gn on time is a symptom of the violation of 
the SEP, but not necessarily of the EEP, whereas the non-constancy of the other constants, 
like the electroweak or strong coupling constants, necessarily represents a violation of the 
equivalence principle in both its forms. 

In this paper, we discuss the effects of time variations of Gn on the light element 
production in Big Bang Nucleosynthesis (BBN), completing and extending the results of 
previous analysis on the subject (see e.g. [HI 111 El) BBN is evidently a good probe of a 
possible time variation of Gn, since it is the earliest event in the history of the universe 
for which we can obtain solid and well-testable predictions. Even a weak (or very peculiar) 
time dependence, which gives no observable effects in high accuracy experiments performed 
at the present epoch, could give sizable effects when translated over cosmological time 
scales. BBN, however, is a complex phenomenon, since each element responds in its own 
way to a modification of Gn- We devote particular attention to this point, introducing 
suitable response functions which relate the abundance of each element to an arbitrary 
time-dependent modification of the early universe expansion rate. 

The plan of the paper is the following: In the next section, we discuss the role of Gn 
in BBN and we derive, analytically, the dependence of the primordial abundances on the 
value of the Newton constant at the key epochs. In section 3, we calculate numerically 
the response functions, emphasizing that different elements are sensitive to the value of the 
Newton constant at slightly different times. In section 4, we discuss the bounds on Gn 
variations that can be obtained from the comparison of theoretical predictions with the 
observational data for light elements primordial abundances. We summarize our results in 
section 5. 

2 The role of in BBN 

The production of light elements (namely ^H, ^He, ^He and ^Li) in BBN is the result of the 
efficiency of weak reactions {p + e ^ n + and related processes) and nuclear reactions 
(which build light nuclei from neutrons and protons) in the expanding universe. The value 
of the gravitational constant determines the expansion rate of the universe and thus, in 

^In this context, the term ’’fundamental parameters” is more appropriate than ’’fundamental constants”. 
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turn, the relevant time scales for the above processes. As a consequence, if we assume that 
the gravitational constant at time of BBN is different from its present value, this translates 
into a variation of light element abundances with respect to standard BBN predictions. 

The above argument clarifies in simple terms the role of Gn in BBN. In order to be more 
quantitative, one needs, as a first step, to identify the key epochs for light element production 
with respect to possible Gn variations. As we shall see, the relevant periods are those 
during which the weak reaction rates and/or the nuclear reaction rates are not vanishing 
nor exceeding the universe expansion rate. Essentially they are the weak-interaction “freeze- 
out” epoch (which occurs at the temperature Tf ~ 0.8 MeV) and the “deuterium bottleneck” 
epoch (which corresponds to Tj ~ 0 . 08 MeV) [ilHlClIin]. In the following, we use GN,f 
and GN,d to indicate the value of the Newton constant during these periods, while we use 
Gn,o to indicate the present value. 

In order to estimate the dependence of the various elemental abundances on Gn^ and 
GN,d) it is necessary to quickly review the basic physical mechanisms responsible for light 
element production. When T ^ Tf, the rate of weak processes which interchange neutrons 
and protons, Tw , is large with respect to the expansion rate of the universe: ^ 

H = , (1) 

where 5 * counts the total number of relativistic degrees of freedom^. As a consequence, 
neutrons and protons are in chemical equilibrium and the neutron abundance = n^jriB, 
defined as the ratio of neutron to baryon densities, is simply given by Xn{T) = [1 + 
exp(Am/T)]“^, where Am ~ 1.29 MeV is the neutron-proton mass difference. 

For T <Tf the weak reaction rate drops below the Hubble expansion rate, the neutron 
abundance freezes out at the equilibrium value X^{Tf) and it then evolves only due to the 
neutron decay: Xy,{t) ~ Xn(Tf) exp(—t/r), where r = 885.7 s is the neutron lifetime. The 
“freeze-out” temperature is basically determined by the condition rw(Tf)/H(Tf) ~ 1 and, 
clearly, depends on the value of the gravitational constant. One obtains: ^ 

/G 

Tf = 0.784 MeV , (2) 

\Gn,o / 

where GN,f is the value of the Newton constant during the freeze-out, i.e. when the tem¬ 
perature of the universe is 0.2 < T < 2 MeV (see next section). The larger is GN,f, the 
earlier is the freeze-out of the neutron abundance, at an higher value, and the larger is the 
^He abundance produced in BBN. 

Light element production in BBN occurs through a sequence of two body reactions, such 
as p(n, 7 )^H, ^H((i, n)^He, ^H((i,p)^H, ^He(p, 7 )‘^He, etc. Deuterium has to be produced in 
appreciable quantity before the other reactions can proceed at all. Nucleosynthesis effec¬ 
tively begins when the rate of deuterium processing through ^H((i, n)^He and ^H(d,p)^H 
reactions becomes comparable with the expansion rate of the universe. By imposing this 
condition one finds the “deuterium bottleneck” temperature which, in standard BBN, is 
given by: Td = 0.08(1 -|- 0.16log(r//10“^°)) MeV 0. 

^Here and in the following, we use a natural system of units in which h = c = ks = 1- 

^Strictly speaking, equation 0 is derived in the context of General Relativity in which Gn is constant. 
In any extension of the standard theory, one has additional terms related to time derivatives of Gn. In this 
paper, we make the usual assumption that Gn is slowly varying (with respect to the early universe expansion 
rate) which implies that these extra terms are negligible. 

'‘See l^rni for the precise numerical calculation of the total weak rate Fw. 
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After nucleosynthesis has started, light nuclei (^H, ^He, ^He and ^Li) are quickly pro¬ 
duced. The ^He abundance is basically determined by the total number of neutrons that 
survive till the onset of nucleosynthesis, since nearly all available neutrons are finally cap¬ 
tured in '^He nuclei. The synthesized elemental abundances of ^H, ^He and ^Li are, instead, 
the result of the complex interplay of the various nuclear reactions efficient during and after 
the d-bottleneck. The value of the Newton constant during this period, GN,d) clearly plays 
a relevant role both for ^He and other elements production. 

The situation with '^He is particularly simple. The primordial ^He mass fraction is 
approximatively given by I 4 ~ 2 Xn(td) — 2 Xn(Tf) exp(—td/r), where is the age of the 
universe at the d-bottleneck. Neglecting the weak dependence of the temperature Tj on 
Gn) one has: 

td ~ 206 (1 - 0.32 log(r?/10-iO)) sec, (3) 

where GN,d is the value of the Newton constant when the temperature of the universe is 
0.02 <T< 0.2 Mev (see next section), i.e. when ^He and the other elements are effectively 
synthesized. By using this formula and considering eq. m, one is able to estimate the 
dependence of I 4 from GN,f, GN,d and rj. One obtains; 

AYa 

dYi = —— = 0.23 5GN,f + 0.09 5GN,d + 0-07 log(? 7 /r 7 cMB) (4) 


where 5GN,i (with i = f, d) represents the fractional variation of the Newton constant at a 
given epoch with respect to its present value: 


<5GN,i 


Gn,! — Gn,0 
Gn.o 


(5) 


and log{r]/ tjcmb) is the logarithmic variation of the baryon to photon ratio with respect to 
the value ?7cmb = 6.14 • 10“^® presently favored by cosmic microwave background (CMB) 
m and deuterium data m- 

The situation with ^H, ^He and ^Li is slightly more complicated. In principle, one has 
to integrate the rate equations, which can be written formally as: 

dy. 

ocrjn^ '^Y xY x {av)T , ( 6 ) 


where Yi indicate the abundance of a given element, the sum runs over the relevant source 
(-|-) and sink (—) terms, and are the thermally-averaged reaction rates.® In order to 

estimate the role of Gn, one can simply note that, since the temperature of the universe 
evolves as dT/dt oc —T^y/G^, the above equation can be rewritten as: 


dYi 7] Hry 

dT “TtG ^ 

'-"N.d 


^ y X y X {av)T , 
-I-,- 


(7) 


®In writing eqs. and Q, we implicitly assumed that only two-baryon reactions are relevant (neglecting 
reactions such asp-te—>n-|-^ and related processes or d -|- 7 —> n -|- p). This is reasonable after deuterium 
bottleneck and provides a good framework to discuss the abundances of ^H, ®He and ^Li, whose abundance 
is essentially established after the d-bottleneck. On the contrary, this is clearly not adeguate to discuss the 
synthesis of '*He. 
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which shows that ^He and ^Li depend on rj and Gn essentially through the combination 
^/^Nd- This suggests that the synthesized elemental abundances of D/H, ^He/H, and 
^Li/H (indicated, in the following, with Y 2 , T 3 , and I 7 respectively) for an arbitrary value 
of GN,d can be related to the standard case (GN,d = Gn,o) through an appropriate rescaling 
in T]-. 

Ti(r/, GN,d) = Tj (^7/(GN,o/GN,d)^^^, Gn,o) • ( 8 ) 

The previous equation, linearized, implies that: 

Ay 

= 7 j(r?o,GN,o) (Alog( 7 ?) - 0.22 5GN,d) (9) 

T i 

where ijq is an arbitrary pivot point, 7 * = (1 /Yi) dYi/d log rj and we have assumed that 
Alog( 7 ) = log(r// 7 o) and <5GN,d are small. We remark that, according to eqs. (I8I9I) . fixed 
values of the abundances correspond in the plane (logr/, (5GN,d) to straight parallel lines 
(with slope ~ 1 / 0 . 22 ). 

The above argument is accurate enough to describe and ^He abundances. However, 
it can be slightly improved in order to predict correctly the behavior of ^Li. In order to 
produce ^Li, one has to use ^He as a target, whose abundance is strongly dependent on 
the value of the Newton constant. This clearly introduce an extra dependence on GN,f and 
GN,d- One expects an extra factor in eq. ((SI) proportional to Y^ which, linearized, gives: 

5Yi = (0.23 (5GN,f + 0.09 5GN,d) +77(70) (Alog(? 7 ) — 0.22(5GN,d) (10) 

where the first terms in the r.h.s are obtained from eq. and we neglected the weak 
dependence of T 4 on the baryon to photon ratio rj. 

The analytical results discussed above allow us to understand the relevance of the various 
physical mechanisms in light elements production and to obtain simple quantitative relations 
(eq. ® for ^He, eq. for and ^He and eq. m for ^Li) between the various elemental 
abundances and the parameters rj, and (^GN,d- In order to check their validity, we 

compare them with the results of numerical calculations. As a first step, we consider 
the case of a constant variation of Gn during the entire period relevant for BBN (i.e. 
<^GN,f = 5GN,d = <^Gn). In this assumption, a linear fit to the numerical result gives for the 
^He abundance: 

5Y^ ~ 0.35 (5Gn + 0.09 log(r//?7cMB) (11) 

with an accuracy at the level of 2% or better in the range (5 Gn = 0.75 — 1.25 and 7 = 

2 10“^® — 10“®, in reasonable agreement with estimate ®. For and ^He , expanding 

around 70 = ^CMB, one obtains: 

5Y2 ^ 72(7cmb) (log(7/7CMB) — 0.25<5 Gn) ( 12 ) 

(5+3 ~ 73(7cmb) (log(7/7cMB) — 0.24(5 Gn) (13) 

in good agreement with predictions ®, with 72 ( 7 cmb) = -3.7 and 73 ( 7 cmb) = -1-3. In 
addition, the ^Li behavior can be described by: 

(5+7 ~ 0.32 (5Gn +77(7cmb) (log(7/7CMB) — 0.22(5 Gn) (I4) 

as predicted by eq. (Uni), with 77 ( 7 cmb) = 4.8.® 

®We remark that, for ®He and ^Li, expanding around an arbitrary value 770 in the range rjo = 

3 • 10 “^® — 10 “® (and using the proper values 77(70)) one obtains the same results as those described by 
relations (HU, (HU and HU, with essentially the same numerical coefficients. 
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3 Response functions 


As underlined in the previous section, the production of each element responds in its own 
way to a variation of the Newton constant. For example, a change of Gn at the 
time of weak interaction freeze-out would have important consequences on the observed 
helium abundance, giving instead negligible corrections to that of deuterium. So far we 
have implicitly assumed that Gn(F) stays constant (at a value not necessarily equal to the 
present one) during BBN, or that it takes two different values at the two key epochs, marked 
as Tf and T^. A more general analysis, that can account for a time dependence of Gn(F') 
along the all BBN period, requires the introduction of suitable functions, which describe 
the response of each elemental abundance to an arbitrary time-dependent modification of 
the early universe expansion rate. 

We have determined numerically the response functions: Qi{r],T), which are defined by: 

6Y,{rj,6H{T)) = 2 J Q,{rj,T) 6H{T) ^ , (15) 


where i = 2,3,4 and 7 and 6H(T) is the fractional variation (assumed to be small) of the 
expansion rate of the universe at the temperature T with respect to its standard value. We 
remark that, in the assumption that Gn(F) is slowly varying, one has that (5 Gn = 2 6H(T) 
from eq. ©• The above equation can then be simply rewritten as 

r dT 

5Y,{r,,5G^{T)) = J g,irj,T) 5G^iT) — , (16) 


which shows that gi{r],T) is basically the functional derivative of lnl 7 (? 7 , (5Gn(T')) with 
respect to 5G^{T). 

The response functions allows us to identify unambiguously the key epochs for the 
production of the various elements and to emphasize that different elements are sensitive 
to the value of the Newton constant at slightly different times. Onr results, calculated 
for r] = rjcMB, are shown in Fig. 1. As expected, the functions gi{r],T) have two peaks 
corresponding to the weak interaction freeze-out and to the epoch, just after the deuterium 
bottleneck, during which the various elements are effectively synthesized. The width of 
the two peaks reflects the fact that the above processes are not instantaneous. One can 
essentially identify the range ATf = (0.2 —2)MeV with the weak interaction freeze-out epoch 
and ATd = (0.02 — 0.2)MeV with the various elements synthesis period. The behavior of 
the functions gi{rj,T) also allows to give a more qnantitative meaning to the parameters 
GN,f and GN,d) which have to be intended, evidently, as the average values of the Newton 
constant during the periods ATf and ATd, respectively. 

The total area under the curves in fig. 1 (integrated in InT) gives the numerical co¬ 
efficient (5Yi/(5GN, which are obtained in eas. fllllllll) in the assumption of constant Gn 
variations. It is interesting, however, to consider separately the early time and the late time 
behavior of the functions gi{r],T) in order to have a feeling of the relative importance of 
the different epochs in the various elements production. In this spirit, we have calculated 
the numerical values: 


Oil 


Pi 


/■2MeV r!T 

/ Q^{r1,T) — 

/0.2MeV 4 

/■0.2MeV UT 

/ Q^{V,T) — 

/0.02MeV J 


(17) 

(18) 
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Figure 1 : The response functions gi{r],T) as a function of the temperature T, for 77 = rycMB- 
The functions describe the effect of an arbitrary time-dependent modihcation of the early 

universe expansion rate on the various elemental abundances (see text for details). 


for i = 2,3,4 and 7, which have to be compared with the coefficients 8Yi/5G-^^i and 
5Yi/5G^^A estimated in eas. (14191101) . For rj = r/cMB, one obtains 02 = 0.12 and [82 = 0.80 
for deuterium, 03 = 0.04 and /Is = 0.29 for hehum-3, 04 = 0.22 and (5^ = 0.12 for helium-4 
and ay = 0.14 and /ly = —0.83 for lithium-7, in good agreement with our predicted values. 

One sees that ^H, ^He and ^Li abundances are essentially determined by the expansion 
rate of the universe during and after the d-bottleneck. As a consequence, the bounds ob¬ 
tained from these elements are basically bounds on (iGN,d- Helium-4, instead, is mainly 
sensitive to the value of the Newton constant during the weak interaction freeze-out. How¬ 
ever, the “response” function Q 4 {T,r]) is rather broad, showing that ^He is sensitive to a 
rather long period, 0.05 < T < 2 MeV, of the early universe evolution. In terms, of 
and GN,d) the bounds obtained from ^He observational data can be considered limits on the 
combination 0.65(5GN,f + 0.35<5GN,d- 

Finally, we remark that the response functions Qi{r],T) have a quite general meaning 
and may be easily applied to discuss any non-standard schemes (e.g. new light - stable or 
decaying - particles, non vanishing muon or tau neutrino chemical potentials, etc.) whose 
main effect is to modify the early universe expansion rate. 

4 The BBN bound on SG]si 

By comparing theoretical predictions with observational data for light element primordial 
abundances one is able, in principle, to obtain bounds for and 5GN,d- However, 

comparison of theoretical results with observational data is not straightforward because the 
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data are subject to poorly known systematic errors and evolutionary effects (see na for a 
review). The present situation can be summarized as it follows: 

i) Recent determinations of deuterium in quasar absorption line systems (QAS) report 
values of D/H in the range D/H ~ 2 — 4 x 10“®. However, the dispersion among the different 
determinations is not consistent with errors in the single measurements. We will use, in the 
following, the value D/H = 2.78lQ3g X 10 ^ given in m, which is the weighted mean of 
most recent deuterium determinations (see m for detailed discussion and references). 
a) Independent determinations of ^He primordial abundance have statistical errors at the 
level of 1 — 2% but differ among each others by about ~ 5%. In particular, by using 
independent data sets, Olive et al. P31ITH] have obtained Y 4 = 0.234 ± 0.003, while Izotov 
et al. jT3 IE] have found I 4 = 0.244 ± 0.002. We will use the “average” value Y 4 = 0.238, 
quoted in m, with the error estimate /YY 4 = 0.005, which is obtained from the dispersion 
of the various I 4 determinations. 

in) The ^Li and ^He primordial abundances are not known, at present, with a level of 
uncertainty comparable to the other elements. The reported values for the primordial 
abundances (see m for ^Li and m for ^He) are, evidently, important as a confirmation 
of the BBN paradigm, but are presently not very effective in constraining possible non¬ 
standard BBN scenarios. For this reason we will not include these elements in our analysis. 

In Fig. 2 we discuss the bounds on (IGn and rj that can be obtained by comparing theoret¬ 
ical predictions with observational data for primordial and ^He. Theoretical calculations 
are made in the assumption of a constant variation, (5 Gn, of the Newton constant in the 
period relevant for BBN. This simple assumption is motivated by the fact that the present 
observational situation does not allow to determine the evolution of GN(t) during BBN. 
However, when considering a theoretical framework in which a specific time-dependence 
for G-^{t) is predicted, one has to keep in mind that ^H and ^He respond differently to a 
non-constant modification of the early universe expansion rate. 

The results shown in Fig. 2 are obtained by defining a x^iV: dG^) as prescribed in m 
which takes into account both observational and theoretical errors in the various elemental 
abundances. The best fit points in the plane (log?],(5 Gn) are obtained by minimizing the 
The three confidence level (C.L.) curves (solid, dashed, and dotted) are defined by 
X^ — Xmin = 2.3, 6.2, II. 8 , corresponding to 68.3%, 95.4% and 99.7% C.L. for two degrees 
of freedom (r/ and (5 Gn), i-e. to the probability intervals designated as I, 2, and 3 standard 
deviation limits. The upper panels are obtained by considering BBN alone, while the 
lower panels show the bounds which can be obtained by combining BBN data with the 
measurement of the baryon to photon ratio from CMB and LSS. This is done by adding 
the contribution: 

X^mb(.) = (19) 

^CMB 

to the BBN chi-square, where ?]cmb = 6.14 • 10“^^ is the baryon-to-photo ratio determined 
by CMB and LSS data and ctcmb = 0.25 • 10“^*^ is the error in this determination j^. In 
this case, the C.L. curves are determined by the condition X^~Xmm — 0, corresponding 

to I, 2, and 3 standard deviations for one degree of freedom ((5 Gn), since g is considered as 
a measured quantity 

The results displayed in fig. 2 allow to obtain the following conclusions: 

^Technically speaking, one should define Xtotfi^GN) = min,,[x^( 7 ], ^Gn) + Xcmb{v)] which depends only 
on the parameter SGn- For graphical reasons and to facilitate the comparison with the bounds obtained 
only from BBN, we presented the results in the two dimensional plane {g, SGn)- 
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i) The and ^He data (panels a and h) select bands in the plane (log(r/), (^Gn), which can 
be easily interpreted in terms of the analitycal (eqs.(jlj) and ©) and numerical leas. ([TT|) and 
d) relations discussed in the previous section. In order to have a bound on (5 Gn one has 
to combine the and '^He observational results (panel c) (in the assumption that 5G^{t) 
stays nearly constant during BBN) and/or to consider the independent information on r/ 
given by CMB+LSS (lower panels). 

a) If we combine and CMB+LSS data (panel d), we obtain (5 Gn = O.OD^Q^fg, in agree¬ 
ment with [ 3 . The quoted bound is consistent with the standard assumption that Gn has 
not varied during the evolution of the universe and is, at present stage, the most robust 
piece of information on the value of the gravitational constant in the early universe. We 
recall that this bound essentially applies to the value of the Newton constant during and 
after the d-bottleneck (i.e. when 0.02MeV < T < 0.2MeV) or, equivalently (in the previous 
sections notations), to the parameter (5GN,d- 

Hi) If we combine the '^He and CMB+LSS data (panel e), we obtain (IGn = —0.11 ± 0.05, 
which shows that '^He observational data favor a reduction of Gn in the early universe, 
even if errors are large enough to allow for the standard value (5 Gn = 0. This result is 
emphasized (reduced) if we consider the “low” (“high”) helium value +4 = 0.234 ± 0.003 
given in 113 GSl (T 4 = 0.244 ± 0.002 given in ^ZHH])) which results in JGn = —0.15 ±0.03 
(JGn = —0.05 ± 0.02) and is also obtained in panels c) and f) where observational 
data are also included. If confirmed (and strengthened) by future data, this could be an 
important indication of non-standard effects in BBN. We remark, however, that the present 
situation is quite delicate. The uncertainty in the quoted bound, (5 Gn = —0.11 ± 0.05, is 
completely dominated by (not well known) systematic errors in ^He measurements. It is 
thus advisable to wait for a better comprehension of these errors, before a final result can 
be obtained. 

iv) In panel /we show the bound, (5 Gn = —0.09 ± 0.05, that is obtained by considering ^H, 
^He and CMB+LSS data. We see that the fit is dominated by helium-4 and CMB+LSS 
observational data, and that the information provided by deuterium only marginally re¬ 
duces the error bar with respect to the previous case. We remark that the above result 
is obtained in the assumption that (5 Gn is nearly constant during BBN. In principle, one 
could fit the data by considering (5GN,f and <5GN,d as independent parameters. In a future 
perspective, this could be interesting as a possible test for variations of Gn during BBN (or 
equivalent schemes). In the present experimental situation this appears too ambitious. The 
value of Xmin — 1-0 foi' the best-fit point indicate, in fact, that the quality of the fit is good 
and that there is no real evidence, at present, in favor of theoretical schemes which predict 
non-constant modification of the early universe expansion rate. 

5 Summary and conclusions 

We summarize the main points of this letter and provide some perspective: 

i) We have discussed the role of Gn in BBN and we have derived analytically the depen¬ 
dence of light element primordial abundances on the values (GN,f and GN,d) of the Newton 
constant at the key epochs: the weak interaction freeze-out epoch and the epoch, just after 
the deuterium bottleneck, during which the elements are effectively synthesized. 

ii) The production of each element responds in its own way to a variation (5 Gn of the Newton 
constant. A general study, that can account for a time dependence of Gn(T) during the 
BBN period, requires the introduction of suitable functions which describe the response of 
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Figure 2: The bounds on SGn which can be obtained from and ^He observational data. The 
upper panels are obtained considering BBN alone. The lower panels are obtained by combining BBN 
data with the measurement of the baryon to photon ratio from CMB and LSS (see text for details). 


each element to an arbitrary time-dependent modification of the early universe expansion 
rate. We have numerically calculated these response functions, obtaining a good agreement 
with the analytical estimates. 

in) We have emphasized that different elements are sensitive to the value of the Newton 
constant at slightly different times. The ^H, ^He and ^Li abundances are essentially deter¬ 
mined by the expansion rate of the universe close to the d-bottleneck. Helium-4 is, instead, 
mainly sensitive to the value of the Newton constant during the weak interaction freeze- 
out. In a future perspective, this could be interesting as a possible test for variations of Gn 
during BBN (or equivalent schemes). 

in) We have discussed the observational bounds on the possible variations of the gravita¬ 
tional constant in the early universe. Our best limit, (5 Gn = O.OOIqjq , is obtained by 
combining observational results with the measurements of the baryon to photon ratio 
obtained from CMB and LSS data. This limit refers to the value of Gn when the tem¬ 
perature of the universe is 0.02 < T < 0.2 MeV (i.e. during and immediately after the 
d-bottleneck epoch) and is consistent with the standard assumption that Gn has not varied 
during the evolution of the universe. 

iv) To conclude, we remark that the results obtained for (5 Gn may be easily applied to a 
more general context in which other constants are allowed to vary and/or new light particles 
are included. The limits on (5GN,f and (5GN,d are, indeed, essentially obtained by compar¬ 
ing the expansion rate of the Universe with the weak reaction rate (in the case of dGyfy) 
or with the light element nuclear reaction rates (in the case of (5GN,d)- One immediately 
understands, then, that the bound on GN,f is basically a bound on ^^(Tf) GN,f/Gp, where 
Gp is the Fermi constant and ^^(rf) is the number of relativistic degrees of freedom at the 
weak interaction freeze-out, while the bound on GN,d is essentially a bound on 5 *(Td) GN,d, 
where g'*(Td) is the number of relativistic degrees of freedom at the d-bottleneck. 
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